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Recent experiments on mesoscopic normal-metal-superconductor heterostructures 
resolve properties on length scales and at low temperatures such that the tempera- 
ture is below the Thouless energy k-&T < Et\i- We describe the properties of these 
systems within the framework of quasiclassical many-body techniques. Diffusive 
and ballistic systems are covered, both in equilibrium and nonequilibrium situa- 
tions. Thereby we demonstrate the common physical basis of various subtopics. 
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1. Introduction 

Quasiclassical techniques have become a widely used tool in the description of mesoscopic 
superconductivity. In comparison to early work which concentrated mostly on the regime near 
the superconducting transition temperature, recent work resolved much smaller length scales and 
covered lower temperatures. At temperatures below the Thouless energy, i.e. in dirty metals of size 
d for k^T < i?xh ~ hD/d 2 , novel features have been observed. An important physical property is 
the proximity effect in a normal metal in contact with a superconductor. It is intimately related to 
the Andreev reflection j^] , a process in which an electron incident from a normal metal tries to enter 
a superconductor. If its energy lies below the gap, it can do so only if it finds a partner electron with 
opposite momentum and spin to form a Cooper pair, leaving a retroreflected hole in the normal 
metal. In this process two charges are transferred into the superconductor. The reflection is phase- 
coherent since the electron and the reflected hole combine in a way consistent with the phase of the 
superconductor. The pair can maintain its phase information in the normal metal over a distance 
~ hvp/2E, where 2E is the energy difference of the electron and hole. This space dependence, 
with diverging length scale at low energies, is responsible for several of the mesoscopic effects to be 
discussed below. 

The proximity effect is characterized by the existence of superconducting correlations in the 
normal metal, expressed by a non-vanishing pair amplitude, F^O. This is to be distinguished from 
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the pair potential A, which in the absence of attractive interactions vanishes in the normal metal. 
Finally, the question whether the normal metal shows a gap in the excitation spectrum is a more 
subtle issue. We will analyze below how the proximity effect reduces the tunneling density of states 
(DOS) at low energies, and, depending on the geometry, induces a gap or a pseudo-gap. 

We will present here a unified view of both equilibrium and nonequilibrium properties of normal- 
metal-superconductor proximity systems. The basis are the quasiclassical Green's functions, which 
will be reviewed in Section 2. Both the Matsubara and the Keldysh formalism are introduced. 
We discuss the clean and dirty limits, which allow systematic approximations. We finish the section 
with remarks on the boundary conditions and the strategies of solution. There exist several excellent 
introductions to this topic [f|-[§ ■ They have been written mostly with the goal of describing massive 
superconductors at temperatures near the transition temperature. Recent experimental work covered 
much lower temperatures and resolved properties on sub-micron length scales. We, therefore, present 
the theory and applications focusing on nano-structured systems at low temperatures. 

In the following sections we will discuss several applications of the formalism and compare the 
results with recent experiments. In the equilibrium case (Section 3) we study the change in the local 
density of states (LDOS) in the vicinity of a normal metal - superconductor interface. It shows a 
strong suppression around the Fermi energy, leading to a mini-gap if the normal metal is finite in 
size 0|. Our results agree with experiments performed by the Saclay group 0. Next, we analyze 
the proximity induced magnetic screening in a normal metal in contact with a superconductor. This 
allows us to describe the experiments of Mota and coworkers on the diamagnetic properties 

of normal - superconductor heterostructures. For a comparison it is necessary to study the induced 
superconductivity in a material with intermediate strength of the impurity scattering. We found 
several new phenomena regarding the non- locality of the current response in this regime [ fL2| . Finally 
we consider the supercurrent through a normal metal sandwiched between two superconductors, 
which allows us to interpret experimental results of Courtois et al. [[l3| . 

In Section 4 we study mesoscopic superconductors and heterostructures in nonequilibrium sit- 
uations. The quasiparticle distribution function can be disturbed, e.g., by external applied gate 
voltages. Under mesoscopic conditions, it acquires a specific double-step form, which has been ob- 
served by the Saclay group fl4| . This nonequilibrium form can be used, in a 4- terminal geometry, 
to tune a supercurrent and, hence, to build a mesoscopic SNS transistor jlB). Moreover, the SNS- 
junction can be switched by the gate voltages to become a 7r-junction. We can explain in this way 
recent experiments by the Groningen group 1 16 - 18 1 . A different SN-heterostructure allows us to 
study the behavior of induced superconducting correlations on a dissipative current. In systems with 
transparent SN-contacts, the dissipative conductance is enhanced, with a peculiar non-monotonic 
temperature- and voltage-dependence, usually referred to as "reentrant" conductance |l9|,|2(|. This 
means, at high temperatures the conductance grows with decreasing T since the proximity effect 
develops, it reaches a a maximum at k^T ss where £/rh is the Thouless energy, and finally it 
drops back to the normal state value at T = 0. This behavior has been observed in experiments 
of the Grenoble group |^l|,|2^]. In more complex, multiply-connected structures the proximity con- 
ductance can be modulated by an external magnetic flux. This modulation is dominated by states 
close to the Fermi surface, for which the spatial range of the proximity effect is large even at fairly 
high temperatures. This makes it possible to detect the influence of the proximity effect over dis- 
tances much larger than the temperature dependent coherence length |l^,|2^,^3|. In systems with 
tunneling barriers, the excess conductance competes with the conductance suppression at the tun- 
neling interfaces governed by the opening of an induced gap. But even without tunneling barriers, 
the conductance may be reduced as a consequence of a 4-terminal measurement in the thin film 
geometry ]24|-|2f|. 

The purpose of this paper is to demonstrate by the use of quasiclassical many-body techniques the 
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common physical basis of several striking phenomena observed in recent experiments on mesoscopic 
superconductors and heterostructures. Of course, neither our presentation nor the bibliography 
can provide a complete review of the field. For a number of further applications of quasiclassical 
techniques we refer to articles in Refs. j2?j-|3^] as well as the articles of Nazarov, Volkov, and Yip in 
this volume. 



2. Quasiclassical formalism 

2.1. Green's functions 

Quantum field theoretical methods in terms of Green's functions are a powerful tool in all many-body 
problems (see e.g. |33|-|35|). Various systematic approximation methods and computation schemes 
have been formulated for them. In this section we will outline, how they can be used, within a qua- 
siclassical approximation to describe mesoscopic normal metal - superconductor proximity systems. 

The starting point for all problems in superconductivity is the Green's function in Nambu space, 
which combines the particle and hole space |56|. Using the pseudo-spinors as a compact notation, 

, we can express the time-ordered Green's function as 



(1) 



The Green's function contains an 'anomalous' component, the pair amplitude F = — i{T^^i), 
characteristic for superconducting systems. Using the BCS Hamiltonian we can write the Gorkov 
equation of motion j33| for G as [J] 

(Go 1 - A-£ imp ) (1,2)® (5(2,1') = 5(1,1') 



(5(1,2)® ((V- A- fW) (2,1') = *(M')- ( 2 ) 



Here 1 and 1' represent sets of space and time coordinates, and ® includes a convolution over the 
coordinates. The free Green's function in Nambu space reads 



G - 1 (l,l') = 6(1-1') 



f3 lr + i^- e0(1)+/i 



(3) 



where dr = Vr — ieA{r)f^ is the gauge-invariant spatial derivative. The part of the electron-phonon 
self-energy which is responsible for superconductivity defines the pair potential 

A(1,1')=J(1,1')(£ ^ with A(1) = A 2 Hm + Ft (2,1). (4) 

Here A is the strength of the attractive interaction. For a normal metal we have A = and hence 
A = 0, but still the proximity effect manifests itself in a non- vanishing pair amplitude F =/= 0. 

To simplify notations, we will from now on use center-of-mass coordinates in space and time and 
Fourier-transform with respect to the relative coordinates. 

Elastic impurity-scattering will be included within the framework of the Born approximation. The 
impurity self-energy then reads 

±(\p\,R,E,t) imp = ^-( I dSG(p,R,E,t)\ , (5) 



2r 



Pr 



f Here and in the following we use units in which Ti = ks = 1. 
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where r is the elastic scattering time, £ = p 2 /(2m) — fi and (. . .) PF denotes averaging over the Fermi 
surface. In contrast, we can neglect inelastic scattering, since we will study mesoscopic samples 
smaller than the inelastic length. Spin-flip scattering can be taken into account by an expression 
similar to (g), see [|-|. 

The Gorkov equation can in principle be used to study mesoscopic proximity systems. However, 
these are genuine inhomogeneous systems. Hence, dealing with full double-coordinate Green's func- 
tions - while not impossible (see e.g. [^7[Q) - may become very cumbersome, and typically one 
makes use of quasiclassical approximations in the course of the calculations anyhow. For a system- 
atic and more efficient approach it is an advantage to perform this approximation already on the 
level of the equations of motion. 



2.2. The quasiclassical approximation 

We will now formulate the equations of motion within the quasiclassical approximation. For sim- 
plicity, with the applications to be discussed in mind, we restrict ourselves to equilibrium and 
stationary nonequilibrium situations, although the quasiclassical approximation can be applied also 
to time-dependent problems, e.g. relaxation processes and collective modes p9[ . 

The Green's function ([!]) oscillates as a function of the relative coordinate \r — r'\ on a scale of 
the Fermi wavelength Ap. This is much shorter than the characteristic length scales in the typical 
problems in superconductivity, £o = wf/A and £t = v-f/2ttT. Moreover, in those problems, it is 
important to study the phase of the two-electron wave function, which according to the definition 
(Q) depends on the center of mass coordinate. For these reasons it is possible, and sufficient for most 
applications, to integrate out the dependence on the relative coordinate. This has been recognized 
first by Eilenberger and by Larkin and Ovchinikov [Q. We can add that the reduction is not 
allowed for another class of mesoscopic effects, e.g. weak localization and persistent currents, which 
are controlled by the phase-coherence of the single- electron wave function, contained in the relative 
coordinates of ([!]). But these effects are usually much weaker than those related to superconductivity. 
On the other hand, the reduction is possible also for problems involving Andreev reflection, since 
the essential information is again contained in the difference of electron and hole wave- vectors close 
to the Fermi surface. 

2.2.1. Gradient expansion 

When integrating over the difference variables the convolution ® in the Gorkov equation Eq. (||) 
(integration over internal variables) requires some care. It can be expressed, after the Fourier trans- 
formation, as a Taylor series 

(A ® B) (p, r , E) = exp % - (d£d* - d£d?) A(p, r, E)B(p, r, E) , (6) 

where r refers to center of mass coordinates. In the problems to be discussed we can neglect short- 
range oscillations, hence we expand this expression up to linear order. 

To proceed, we first subtract the Gorkov equation (||) from its conjugated form. We observe, that 
the Green's functions and the self-energies arc linear combinations of Pauli matrices not including 
the unit matrix. This simplifies the equation of motion to 



ipd r + Ef 3 - zA - E, G 



|9 r (A + t + e4> + m), V P G} + |v p S, d r G*} = . (7) 



This form is much simpler than the original. We note that it still accounts for particle-hole asym- 
metry, which is necessary, e.g., for the description of thermoelectric effects p9|,|40|. 
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2.2.2. Quasiclassical Green's functions 



On the one hand, we want to ignore the information contained in the fast oscillations of the full 
Green's function G as a function of \r\ — r%\, which produces after Fourier transformation a pro- 
nounced peak at |p| = pp. On the other hand, we have to pay attention to the dependence on 
the transport direction, i.e. on the direction of the velocity, Vp, at the Fermi surface. To make this 

2 

explicit we write G(£,,vp,r, E), where £ = fc- — fi depends on the magnitude of the momentum. 
The quasiclassical Green's function is then defined by 



if 



g(r,vp,E) = - J-d£G(t,v F ,r,E) , 



(8) 



where the integration contour has two parts covering both half planes, see Q for details. From 
Eq. (Q) we get, from now on setting p = p-p, the Eilenberger equation of motion for quasiclassical 
Green's functions, 



v F d,g(r, v F ,E) 



-iEf 3 + A 



1 

27 



(g(r,v F ,E)) VF ,g(r,vp,E) 



The elements of g are 



(9) 



(10) 



which is still a linear combination of three Pauli matrices fi/2/3 with /t being the time-reversed 
counterpart of /. This symmetry can be used for convenient parameterizations of the Eilenberger 
equation (see Section 2.6). 

As the right-hand side of Eq. (^) vanishes (in contrast to the Gorkov equation (||)), it only defines 
the Green's function up to a multiplicative constant. The constant can be fixed by the following 
argument: As g is a linear combination of the Pauli matrices, the square of the Green's function 
is proportional to the unit matrix gg = cl. From Eq. (||) we see that the proportionality constant 
c is space- independent. If we now consider a system containing a sufficiently large superconductor, 
we can identify a region "deep inside the superconductor" . Here g equals its bulk value, which 
can be calculated by performing the steps used in the quasiclassical approximation for the known 
solution of the homogeneous Gorkov equation Eq. (||). This procedure yields c = 1, i.e. the 
Green's functions are normalized 

gg=l, i.e. g 2 + ff = 1 . (11) 
A more general, mathematical derivation of the normalization condition can be found in p9[ 



2.3. Kinetics and time-ordering 

From the knowledge of the retarded and advanced Green's functions alone, we can calculate energy- 
dependent quantities like the density of states 

N(r,E)=N Q Re{(g R (r,v F ,E))}. (12) 

This, and similar expressions including the off-diagonal Green's function will be denoted in the 
following as spectral quantities. In thermal equilibrium these quantities also determine all properties 
of our system. Equivalently, the thermal Green's function maybe calculated in Matsubara imaginary- 
time technique. Then, the spectral quantities are found by analytic continuation. 
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In general, we need in addition information on how the quasiparticle states are occupied, i.e. about 
distribution functions. To evaluate those under noncquilibrium conditions we will use the Keldysh 
technique. Both techniques are described in numerous references, so we will not rederive them but 
rather explain the basic idea and provide a practical guideline for using them in various situations. 

2.3.1. Matsubara technique 

The Matsubara Green's functions technique |ll| has been developed to describe many-body systems 
in equilibrium at finite temperature fl33|, p4[. In thermal equilibrium, the eigenvalues of physical 
observables do not depend on (real) time. In order to calculate thermal expectation values, we 
trace over all states using the Boltzmann weight exp(— -ff/T), which can be viewed as analytic 
continuation of the time-evolution operator exp(iHt) to imaginary direction r = it. It is sufficient 
to know this operator (and with it the Green's function) only in the interval < r < l/T. By Fourier 
transformation and exploiting the Fermionic symmetry one sees that all the necessary information 
is contained in the Green's functions defined for a discrete set of energies E = iu) n , proportional to 
the Matsubara frequencies uj n = [2n + \)itT with integer values of n. 

For those frequencies, in a bulk superconductor, the Green's functions read 

, 9uJ = uj/Vt , f u = A/fl with n = V|A 2 | +lu 2 . (13) 

This form will serve in the following frequently as a boundary condition. The Green's functions in 
imaginary times show usually no singular or oscillatory structure and behave rather smooth and 
monotonic. This, together with the fact that is real in the absence of fields or phase gradients, 
simplifies the numerical calculations in thermal equilibrium. 

Expectation values of physical quantities can be expressed via Green's functions. To calculate 
thermal averages we rotate onto the imaginary time axis and perform the quasiclassical approxima- 
tion. For instance, the result for the supercurrent density is 

j(r) = -2ie7dV T^Tr(v F f 3 5c(wF,r-)) . (14) 

UJ 

In general, it depends non-locally on the applied fields. Similarly, the self-consistency equation for 
the pair potential can be obtained from the corresponding expression (|j) for Gorkov's Green's 
functions ]33]] . The integrations leading to the quasiclassical Green's function can be performed, and 
the self-consistency relation reads 

A = XN 2ttTY / (D- (15) 

2.3.2. The Keldysh technique 

The Keldysh Green's function technique Q allows describing many-body systems outside equilib- 
rium. Pedagogical reviews are given in articles || and books [ ff3[ . It has proven particularly useful 
in the description of nonequilibrium superconductors ^,|^, ^9|. Again, we will not rederive the for- 
malism, rather we summarize the concepts and some major results. We will later restrict ourselves 
to the dirty limit, although the Keldysh technique can be applied to ballistic systems as well |44| ]. 

By means of the Keldysh technique it is possible to describe the real-time evolution of systems 
in noncquilibrium and at finite temperature. We assume that the system is initially, at t = — oo, in 
a thermal equilibrium state. The time evolution of the system, for instance of its density matrix, is 
described by a forward and a backward propagator. Keldysh showed that this real-time evolution, 
as well as an evolution along the imaginary time axis to account for thermal averaging, can be 
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combined to a propagation along a single contour in the complex time plane, with a forward and 
backward branch running parallel to the real axis and a vertical part between t = — oo — i/T and 
t = — oo J42] ]. Green's functions are defined on this contour with time-ordering along the Keldysh 
contour. Depending on whether the time arguments are on the forward or backward part, the Green's 
functions reduce to different analytic parts. For instance, the time-ordered Green's functions or the 
Kadanoff functions G < ^ > - ) can be obtained 35 1. From those one obtains by appropriate linear 
combinations the retarded and advanced Green's functions G R ' A - ) as well as the Keldysh Green's 
function G K . They are related by 



6>(l,l , ) = *<* t (l , )*(l)) , 6<(l,l , ) = -<(*(l)* t (l / )), 



(16) 



and 



G R (l,l') = 0(ii-ti) G<(1,1')-G > (1,1' 
G A (1,1') = -9(t\- tl )[G< (1,1') -G>(1,1>) 
G K (1,1') = G > (1,1') + G < (1,1'). (17) 
To compactify notations, these Green's functions can be written in matrix form in 'Keldysh space', 



G 



G R G K 
G A 



(18) 



If we describe superconductors, each entry of this matrix is still a 2 x 2 matrix in Nambu space. 

The retarded and advanced Green's functions G R ^ determine what we will call spectral, i.e. 
energy-dependent properties of the system. Their Eilenberger equations have just the standard 
form, viz., Eq. (||). Usually they can be obtained by analytical continuation of the Matsubara 
Green's function if we set u> — ► — iE±0. On the other hand, G K is needed to account for properties 
of the system which depend on the noncquilibrium distribution function. We will later relate it to 
this physical quantity and indicate ways to calculate it. Expectation values of physical quantities 
can again be related to the Green's functions, now defined on the Keldysh contour |35[| . For 
instance the electrical current becomes 



j(r) = -eN Q / dE Tv(v F f 3 g K (v F , r, E)) 



(19) 



2.4. The dirty limit 



2.4.1. Usadel equations in Matsubara technique 

Frequently the superconducting material has strong impurity scattering and is described by the 
so-called 'dirty limit'. The requirement is that the elastic scattering self-energy dominates all other 
terms in the Eilenberger equation. In this limit the electron motion is diffusive and the Green's 
functions are nearly isotropic. I.e. one can expand the Green's functions in spherical harmonics 



9u(r, v F ) = G u (r) + v F g u {r) 



(20) 
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Following the prevailing convention we denote the angular average of the quasiclassical Green's func- 
tion again by a capital letter. In the expansion we assumed that v^g <C G. From the normalization 
condition (pi]) it follows that G^(r) = 1 and {G w (r), g w (r)} = 0. Angular averaging of (0) yields 



wf 3 + A(r), G w (r) 



while averaging of (^|) after multiplication by t>p, yields 

9 w ( r ) = - T Gu>{r) d r , G w (r) 



(21) 



(22) 



Here the condition g(r)/r ujt^ + A has been used. Inserting ( f22| ) into ( |2l| ) leads to the Usadel 
equation [49], 



D 



wfs + A(r), G u {r) 



(23) 



where D = Wpr/3 is the diffusion constant. This equation is much simpler than the original Eilen- 
berger equations and has been widely applied to describe properties of mesoscopic proximity systems. 
The current becomes in the dirty limit 



j{r) = -^ r E Trf 3G.(r) \8 r ,G u (r) 



(24) 



where ctn = 2c 2 NqD is the normal-state conductivity. Note, that a additional restrictions on the 
mean free path are necessary for this local relation to hold. We will discuss this in Section 3.2. E.g. 
for a bulk superconductor with a real order parameter it is given by 



A 2 

j(r) = -wa N T J2 7j2 A (r) , 



(25) 



which demonstrates an important property of the dirty limit, namely a local supercurrent-field re- 
lation. 



2.4.2. The Usadel equation for Keldysh Green's functions 



The reduction to the dirty limit can be performed, similar to the procedure outlined above, also 
for the Keldysh Green's functions. The Usadel equation in Keldysh x Nambu space thus reads 



D 



B r ,G(E,r) d r ,G(E,r) = [-iEf 3 + A, G(E, r)] 



where 



r 3 



r 3 
f 3 



A 



A 
A 



(26) 



(27) 



Although the Keldysh technique works also for time-dependent situations we restrict our- 

selves here to stationary nonequilibrium problems. We further consider only structures in which all 
superconducting reservoirs are at voltage V = 0, thus avoiding effects related to the nonequilib- 
rium Josephson effect. The normalization condition still holds in the form GG = 1. In terms of the 
components of ( jl8f ) it implies 



qRqR = qAqA = l ^ and qKqK + qKqA = q 



(28) 
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As a consequence of the second relation G K can be parameterized as 

G K = G R h - hG A . 



(29) 



From the Kcldysh component (i.e. upper right) of the Keldysh-Usadel equation (26) we obtain a 
kinetic equation of motion for the distribution matrix h 



D 



V 2 h + (G R VG R )Wh - V/i(G A VG A ) - V(G R (Vh)G p 
- (G R [A, h] - [A, h]G A ) + iE (G R [h, f 3 ] - [h, f 3 ]G 



= 



(30) 



Here, we made use of the fact that G R ^ satisfy the respective components of the Usadel equation. 
The kinetic equation has only two independent entries. Since, furthermore, it is a linear equation, 
we can assume h to be diagonal. Returning to the definitions, we can relate it to the distribution 
functions for electrons and holes 



1 - 2fei 
2f h - 1 



(31) 



where the energy is measured from the chemical potential of the superconductor. In thermal equi- 
librium, for instance in the reservoirs at voltage Vr, it can be expressed by Fermi functions 

_/tanh(^) 



1 - 2f(E) 






2f(-E) - 1 





tanh(^§J^ 



(32) 



For practical calculations, it is convenient to split the distribution matrix into odd and even com- 
ponents with respect to the Fermi surface h = /l + /t^. The notation / L / T refers to 'longitudinal' 
and 'transverse' changes of the order parameter associated with the respective nonequilibrium dis- 
tributions 1 46 , 4^] . Here it is sufficient to mention that deviations of /l from equilibrium are related 
to an effective temperature change, whereas /t is related to an effective chemical potential shift. 
The equilibrium forms in the reservoirs are 



/l(T) - 2 



tanh 



E + Vr 
2T 



+ (-) tanh 



E-Vr 
2T 



(33) 



We can recover the usual electron distribution function in the end of the calculation as 2f(E) 
l-f h (E)-f T (E). 

The kinetic equations for the two components reduce to two coupled diffusion equations 



V(X>TV/ T )+2Im{j E }V/ L = 2^/ T 
V(P L V/ L ) + 2Im{j B }V/ T = 
with energy dependent spectral quantities 

f 3 (G R VG r - G A VG A 



2lm{j E } 
K 



iTr 

4 



jTr ( 1 - G R T 3 G A r 3 



-Tr ( 1 - G R G A 
4 



= — Tr 



A (G R + G A ) 



(34) 
(35) 

(36) 
(37) 
(38) 
(39) 



Note, that the diffusion constant is generalized to energy dependent spectral diffusion coefficients, 
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D T / L (i?). They are even functions of energy, whereas the the spectral supercurrent, jg, is an odd 
function of the energy. Recently, it has been argued, that in some particular situations, the diffusion 



equation should be supplemented by another term usually neglected 48 1. 

The right-hand side of eq. ( |34| ) implies, that /x relaxes to inside the superconductor on a scale 
of £o: i- e - it adjusts the chemical potential to that of the superconductor. In contrast the right- 
hand-side of Eq. (^) vanishes. This component /l relaxes only due to inelastic processes, which are 
neglected here, on a longer length scale 0|. 

2.4.3. Physical quantities and Keldysh Green's functions 

We now list some results for physical quantities expressed by Keldysh-Usadel Green's functions. 
The self-consistency equation for the order parameter reads in analogy to ( fi"5| ) 

A = ^y dEf L (F R ~F A ) , (40) 

while the current becomes 

3 = T / dETl [^^VG* + G K VG A )] . (41) 

We can now use the parameterization of the Keldysh Green's function by distribution functions, 
Eq. (^9|), and obtain terms proportional to the two distribution functions and their respective 
derivatives. By studying the parity of the terms in energy we note that some terms vanish, and we 
end up with the supercurrent 

/oo 
JE f L (E)lm{j E } (42) 

and the dissipative current component 

Jn = ^-J dE V T {E) Vf T (E) . (43) 

The kinetic equation ([54]) guarantees the conservation of the total current, j = j s + j n , in the nor- 
mal metal. The other kinetic equation ( p5| ) can be identified as describing the convective part of the 
thermal current. 



2.5. Boundary Conditions 

In the quasiclassical approximation the information on length scales of the order of the Fermi 
wavelength have been integrated out. Consequently one can no longer account directly for the effect 
of potential barriers or interfaces on this level. It turns out, however, from a study of the full theory, 
that potential barriers and interfaces can be accounted for by effective boundary conditions for the 
quasiclassical Green's functions. A derivation of boundary conditions valid at arbitrary transmission 
of the interface has been given by Zaitsev ]49| ]. These boundary conditions couple the classically 
transmitted and reflected trajectories. The corresponding wave vectors k 1 2 are shown in Fig. ^. The 
in- and outgoing Green's functions gf 2 on both sides (1 and 2) on these trajectories are combined 
as indicated in Fig. [l]. They fulfill the boundary conditions 

9t = 92, 

t[R(gt) 2 + (9s) 2 } = Tg-g+, (44) 
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x = 



g+ =g(0~,k+) 



Si = 9(0 ,fci ) 



9i' a = l(9i ±9i ) 




g+=g(0+,fc+) 



9 2 =9(° + > fc 2 ) 



\{9i ±9 2 ) 



|(9 S l±9l) 



Fig. 1. Trajectories at an interface. The parallel momenta are conserved, but for different Fermi velocities on both 
sides the perpendicular momenta differ. The Green's functions on the trajectories involved are indicated by g\2~ ■ 



where R = 1 — T is the momentum-dependent reflection coefficient. The functions gf 2 and g^'~ are 
defined in Fig. [l]. Note, that current conservation through the boundary is ensured by the first of 
these equations. The boundary conditions are also valid for Keldysh-Green's function, provided the 
Nambu matrices are replaced by Keldysh matrices. 

The boundary conditions can be simplified in the dirty limit for low transparencies of the inter- 
faces. As was shown by Kuprianov and Lukichev |50f they read in this limit 



Here t — {py2xT ' / 'pp^R) is a (small) parameter related to the transparency, and py2x is the projec- 
tion of the Fermi momentum onto the normal of the interface. Recently Lambert et al. plj ] showed 
that the second of these conditions constitutes the first term of an expansion in the parameter t, 
and they calculated the second term in this expansion. 

2.6. Parameterizations 

For the further handling of the quasiclassical equations, both numerically and analytically, two 
parameterizations have turned out to be especially useful: the Riccati- and the ^-parameterization. 

2.6.1. Riccati parameterization 

We will make use of the Riccati parameterization mostly in equilibrium problems for general strength 
of the impurity scattering. In this case we write the Matsubara Green's functions as j52| 



ax " ax 

I2G2-J—G2 = t Cr2,Gl 



(45) 




(46) 
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where the functions a w {r,v-p) and (^(r^up) obey the equations 

-v F Va = 2uja + A*a 2 -A, 
v F Va? = 2cDa t + Aa f2 - A* . (47) 

Here 

<D=w+ (g) /2r and A = A + (/) /2r (48) 

are the renormalized energy and pair potential. These equations are of the Riccati type. The func- 
tions a and a) appearing in this parameterization are directly related to the coefficients of the 
Andreev amplitudes. This can be seen, by noting that (^) in the clean limit directly follow from 
the Andreev equations Q identifying a — u/v ]53( ]. For numerical purposes these equations are 
well suited, since they are two uncoupled stable differential equations (one for each direction of in- 
tegration). This should be contrasted to the original Eilenberger equations, which are coupled and 
unstable differential equations. These equations also provide a basis to treat the linear response to 
an external field 

2.6.2. 9 parameterization 

The dirty-limit equations, which do not dependent on vf any more, allow a simpler parameterization 

G R (E r) = ( cosh(6>) sinh(0) exp(ix) 

^ ' ' I — sinh(6>) exp(— i\) — cosh(#) 

G A (E, r) = ( . -™ sh ^ \ . (49) 

v ' ' \ — smh(f) exp(— v\) cosh(0) J y ' 

Here, 9(E, r) and x{E, r) are complex functions the bar denotes complex conjugation Q The Usadel 
equations in the normal metal are now written as 

DdlO = -2i£sinh(9+ y (d x xf smh26 

d x j E = , ^ = 2sinh 2 ^ X . (50) 
In imaginary time, it is more convenient to parameterize |55f 

~ cos(6») sin(0) exp(ix) \ 

9uj[T > ~ \ sin(#)expHx) -coe(fl) )' [ ' 

where 9(uj,r) and x(w,r) are now real functions. At the interface to a superconducting reservoir 
with pair amplitude Age 1 ^ (which vanishes in normal conductors), which is connected to the probe 
by good metallic contacts, the boundary conditions read 

9 = 6 S = Arctanh(A/£) , 
X = 0o ■ 

For a normal reservoir, we set 9 = 0. In the presence of tunneling barriers, the boundary condition 
get modified. For instance, if we consider a normal layer with resistance per length i?N coupled via 
a barrier with resistance Rt to a superconductor, the condition on 9 following from Eq. (^|) reads 

r t dj = sinh(0 - 9b) . (53) 

f It is also possible to parameterize using sin- and cos-functions. We prefer the present convention, which 
requires to choose the gauge of A such that for a real pair potential, A = |A|, we have A oc ti. Compared to 
our previous definition in eq. (0) we have to shift the phase of the pair potential, <f>, by a gauge transformation 
of the order parameter phase Dy tv/2. 
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The direction of the derivative points away from the superconductor described by 6s, and r t = 
Rt/R>N is the ratio of the resistances. For the kinetic equation Eq. (|35|), the boundary condition is 

n Dd x f T/L = M t /l(/to/lo ~ /t/l) . (54) 



where /to/lo are the components of the reservoir distribution function (33). Using the Green's 
functions G5 and Fg 1 of the superconductor we defined the spectral functions 

M T/L = Re(G^)Re(G R ) ± Im(F|')Im(Ft R ) . (55) 



2.7. Solution strategies 

In some cases we can derive approximate analytical solutions, for instance linearize the equation 
(see e.g. pq|). In general, however, we have to rely on numerical solutions. To find the solution for 
a given transport problem, we can proceed as follows: 

1. Start with a given A(r). 

2. Solve the retarded Usadel equation. 

3. Take the solutions to calculate N(E), Je, and T> L / T . 

4. Solve the kinetic equations. 

5. Calculate a new A(r), and iterate until self-consistency is achieved. 

It is worth noting that considerable simplifications arise if one is only interested in equilibrium prop- 
erties. Then it is possible to determine the self-consistent problem in the Matsubara representation, 
which is numerically much more stable than the real-time version. (In the latter one encounters 
frequently the cancelation of large terms.) Even if one is interested in transport properties it is of 
advantage to proceed in two steps. First, the pair potential, which does not depend on energy, is 
determined in the Matsubara representation. Then, after self-consistency has been achieved, one 
proceeds to solve the real-time/ real-energy quasiclassical equations to obtain spectral quantities 
and the transport properties. 



3. Spectral quantities and equilibrium properties 
3.1. Density of states 

As a first application we study of the local density of states in spatially inhomogeneous dirty systems. 
Examples are a normal metal in contact with a superconductor, inhomogeneous superconductors 
or boundaries of superconductors with unconventional symmetry. Depending on the problem, the 
pair potential A(r) has to be determined self-consistently. As mentioned above, for the ease of the 
numerical procedure this is best performed in imaginary frequencies. The next task is to solve the 
Usadel equation for real energies for a given pair potential. Once we have determined the matrix 
Green's function (j(E,r), the local density of states is given by 

N(E, r) = N Re \g r (E, r)l = N Re [cosh(6(E, r)} , (56) 



and similar for other spectral functions. 
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In this chapter we will calculate the local density of states N(E,r) for geometries in which 
a normal metal and a superconductor are in contact [^j . First we consider the simplest case of an 
infinite normal metal in contact with an infinite superconductor. In this system the spatially resolved 
density of states has recently been measured with the help of tunnel junctions contacting the normal 
metal in different positions ^. As a second example, we consider a finite thickness normal metal in 
contact with a superconductor. In this case, there is a gap in the excitation spectrum of the normal 
metal. The form of the density of states in this system leads to rather unusual current-voltage 
characteristics of SNINS tunnel junctions. This has recently been measured in proximity point 
contacts 57 , in which the 'insulator' is an atomic point contact with variable transparency f38f . 



The final system examined is an overlap junction, where the normal metal and the superconductor 
have an finite overlap. Interestingly, in this system a local gap-like feature in the overlap region 
exists. In all examples in this section the pair potential is calculated self-consistently. 

3.1.1. Infinite proximity system 

We consider a quasi one-dimensional geometry, i. e., an infinite normal metal is in contact with a 
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Fig. 2. Density of states on the normal side of an N-S boundary for two spin-flip scattering rates: r s f = (dotted 
lines) and V s { = 0.015A (solid lines). 



superconductor of the same cross section. The solution of the Usadel equations yields the DOS on 
the normal side at different distances from the NS-boundary as shown in Fig. ||. The peak height 
and position change with distance. In the absence of pair-breaking effects the DOS vanishes at the 
Fermi level for all distances (dotted curves). Inclusion of a pair- breaking mechanism, T s f ^ 0, (solid 
curves) regularizes the DOS at the Fermi level, and also suppresses the peak height. The curves are 
in qualitative agreement with experimental data presented in Ref . |^] . The self-consistent calculation 
presented here leads to a slightly better fit than the curves derived in Ref. M, in which a constant 
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pair potential was assumed. In particular, the low-energy behavior of the experimental curves is 
reproduced better. 

3.1.2. Minigap in a finite normal metal 

Next we consider a normal layer with finite thickness ~ £s in contact with a semi-infinite 
superconductor occupying the half space x < 0. This system displays what is called a 'minigap' in 
the density of states of the normal metal. This gap has been first noted by McMillan J5?| within 
a tunneling model which ignores the spatial dependence of the pair amplitude. We consider here 
the opposite limit, assuming perfect transparency of the interface but accounting for the spatial 
dependence of the Green's functions. 



EVA 




Fig. 3. Minigap E g as a function of the normal-layer thickness. Inset: local DOS of an N-layer of thickness 
djq = l.l£s m proximity with an bulk superconductor . 



The boundary conditions at the interface to the superconductor remains unchanged as compared 
to what has been discussed above. But at the open side, x — g?n, the condition is d9(E, x)/dx = 0, 
which is the same as if the normal metal were bounded by an insulator. In this case the DOS in 
the normal layer develops a minigap at the Fermi energy. It is smaller than the superconducting 
gap, and decreases as the thickness of the normal layer is increased. Results obtained from the 
self-consistent treatment are shown in Fig. |^. Details of the shape of the DOS depend on the 
position within in the N-layer jij. However, the magnitude of the minigap is space- independent, as 
displayed by the inset of Fig. ||. The magnitude of the gap is expected to be related to the Thouless 
energy £/rh ~ D/d^, which becomes a relevant energy scale in mesoscopic proximity system. This 
relation has to be modified in the limit c?n — > 0. Indeed, as shown in Fig. pi a relation of the form 
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E g ~ (const £s + ^n) -2 fits numerically quite well. The sum of the lengths may be interpreted as 
an effective thickness of the N-layer since the quasiparticle states penetrate into the superconductor 
to distances of the order of £g . The effect of spin- flip scattering in the normal metal on the minigap 
structure is also shown in the inset of Fig. |3|. The minigap is suppressed as T s f is increased until it 
vanishes at T s f w 0.4A. For T s f — our results for the structure of the DOS agree with previous 
findings A minigap was also found in heterostructures with low transparency barriers [ pT[ 

and in a two-dimensional electron gas in contact to a superconductor |6^] . A minigap in the DOS 
has also been found in a quantum dot coupled to a superconductor, if the shape of the dot is chosen 
such that it shows chaotic behavior in the classical limit fS3| . This systems are frequently described 
by random matrix theory. In contrast, dot structures which are integrable in the classical limit have 
a non-vanishing DOS at low energy. This behavior is what one finds in clean normal metal coupled 
to a superconductor Q . 



3.1.3. Overlap junction 

Finally we consider a geometry, in which the normal metal and the superconductor have a finite 
overlap region in x-direction (see Fig. ^) which is longer than £g. This geometry tries to mimic 
the experimental realization of junctions between different metals. A measurement of the density 
of states along such a junction would give additional information on the details of the proximity 
effect. In principle the problem is two-dimensional. But if the total thickness ds + djy of the two 
layers is less than the coherence length, we can neglect all derivatives in the transverse direction and 
simply average the Usadel equation over the width. This result for the self-consistently determined 
pair potential is shown in the lower graph of Fig. ^. In the middle of the junction the pair potential 
saturates at 0.5 As- 



N(E,x)/N 
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Fig. 4. Geometry and pair potential of an 
overlap junction for ds/(dg + d^f) = 0.7. The 
pair potential saturates around 0.5As in the 
middle of the junction. 




Fig. 5. Density of states along an overlap junction. The curves for different positions are offset by 0.1 on the 
vertical axis. The distance between two adjacent LDOS is 2.5£s and the position is indicated by the bar on the left 
hand side. One clearly sees the BCS-like density of states with reduced gap in the middle of the overlap junction. 
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The density of states along the overlap junction develops rather interesting structures as shown 
in Fig. U The DOS evolves gradually from the BCS-behavior in the lowest curve to the normal 
metal behavior in the topmost curve. In the middle of the junction a BCS-like peak in the density 
of states appears at an intermediate energy due to the aforementioned plateau in the pair potential. 



3.2. Magnetic response of a proximity structure 

The induced superconductivity enables a normal metal to produce magnetic screening [ |65| , |S6[ 
reminiscent of the usual Meissner effect in superconductors. However, it depends on mesoscopic 
scales, like the thermal coherence length £t = vf/2ttT, and the thickness of the normal metal 
g?n fit] lll-f^] 0- In addition, a length related to the penetration depth of superconductors, the 
London lcnghtA^j = 4ire 2 n e /m plays an important role. 

In the following we we develop a linear response relation between the current and the vector poten- 
tial in terms of the zero-field Green's functions in the presence of arbitrary impurity concentration. 
We will show, that this formula reproduces previous results for the magnetic screening obtained 
in the clean and dirty limit. Finally we show that for intermediate strength impurity scattering 
qualitatively different results arise, consistent with the observations in recent experiments. 

The nonlinear response of the normal metal under proximity shows interesting properties too 
f35|, |73t [74| . At low temperatures, the magnetic response exhibits a first order transition both in the 
clean and in the dirty limit. At a certain critical field, called the breakdown field, the induced mag- 
netization is suddenly reduced, leading to a less diamagetic state, or even vanishing diamagnetism 
in the clean case. Both, the breakdown field and the size of the jump depend on temperature. All 
these features have been measured experimentally [ fL0| . The qualitative agreement of the temper- 
ature dependence of the breakdown field with the clean limit theory |Q showed that the samples 
are close to the ballistic regime. 

3.2.1. Linear response formalism 

We consider the magnetic screening in linear response of a system composed of a superconductor 
occupying the half-space x < 0, in contact with a normal layer of finite thickness occupying the 
region < x < g?n, to a magnetic field B(r) — B(x)z in z-direction. We chose a gauge such that the 
vector potential is A(r) = A(x)y. For the analysis it is convenient to expand the Green's functions 
in the Riccati parameterization 

a>u) {v x > v y , x) = <z° (v x , x) + iai (v x , v y , x) (57) 

and a similar expansion for a) . Here a and a 1 are the first terms is in an expansion in A, and 
we have explicitly indicated the dependence on the two relevant components of vp. Accordingly we 
expand the diagonal and off-diagonal Green's functions g u and In first order the Eilenberger 
equation reads 

v d r ~ 1 

- Yd^ a ^ Vx,v y ,X ^ = + A ( x ) a Z (v x ,x) al{v x ,v y ,x) + ev y A(x)a^(v x ,x). (58) 

The symmetry relation a* (x^, v y , x) = — a^(ur, — v y , x), which is a consequence of the gauge chosen, 



implies that the impurity self-energies contained in ui and A, see Eq. (|4J), depend on a° only. Thus, 
we can integrate the linear inhomogeneous relation ( |5^ ) directly. The solution still depends on the 
boundary conditions. Assuming perfect transparent N-S boundaries and specular reflection at the 
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boundary to the vacuum we obtain for the current 



3y{x) = - 



2 9 



TT 



u>0 



dx' I dv 



JL[l + gl(v x , X )][l-gl(v x , X >)] 



<d(x - x')m u {v x ,x, x 1 ) + Q(x' - x)m u (-v x ,x, x') + m u {-v x ,x, du)mu)(v x ,dn, x') 



(59) 



A(x') 



K(x,x')A(x')dx' , 



where g^{v Xl x) is related to a^(v x ,x) by (|46|), and 



mu;(v x ,x,x') 



exp — 



A(x") 
0(v x ,x") 



-dx" 



(60) 



The relation (^9|) is characterized by the interplay between the range of the kernel determined by 
the decay of m u (v x ,x, x') and the amplitude of the prefactor + — g^{x')\. This determines 

the non-locality of the current- vector potential relation. 

To examplify the usefulness of this result we reproduce the current response of a half-infinite 
superconductor. Setting d^ = 0, the solution of the Eilenberger equation (J9j) takes the simple form 



g a = w/fi, f = 4 = A/O, where Q = (A 2 
response kernel 



2 ) 1 / 2 . Inserting this in ( |59j ) we obtain the linear- 



K s (x,x') 



e 2 p F 



jdu 

LU>0 ^ 



l — u 2 /v 



£ [ e -(20+l)J^i_ 



3 (2fi+l)^±2L 



(61) 



This result describes the current response of an arbitrary superconductor, as first derived by Gorkov 
[jj3f , which here additionally includes the effect of the boundary. For fields varying rapidly spatially 
we arrive at a non-local current-field relation of the Pippard-type (75) , while for slowly varying fields 
the kernel can be integrated out in Eq. (|59|) , producing the local London result Q . 

The solution ( |59| ) is valid rather generally. We now concentrate on a structure with a normal 
metal of thickness di$ in contact with a semi-infinite superconductor. For simplicity we neglect the 
penetration of the field into the superconductor by choosing A(x = 0) = 0. This approximation 
leads to corrections to the induced magnetization of the normal metal of order As/c^n only. As the 
second boundary condition we put dA(x) / dx\ x=c i N = H with the applied magnetic field H. The 
spatial integration is now restricted to [0, g?n] and m takes the form 



m u {v x ,x,x') = exp ( — / 
V V xT J x 



fl\v x ,x") 



-dx" 



(62) 



Here, (f°(x)) = j^duf^vpu.x) and f°(v x ,x) is related to a°(v x ,x) by (|||). 

The magnetic response of the normal-metal layer is expressed by the magnetic susceptibility 

A(d N ) 



X 



1 

47T 



Hd N 



(63) 



It's calculation requires a self-consistent solution of Eq. (^9|) and the Maxwell equation, 
d 2 A(x)jdx 2 = —4:Trj y (x). We will now present results for the clean and dirty limits, as well as 
the intermediate regime. 
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3.2.2. Clean limit 

In the clean limit, r — > oo, and consequently m — 1. Furthermore, <?JJ( u :c) is constant in space, and 
the current takes the form [fj7f 



where 



M - = _ 47rA 2 (T) <2 n / A{x)dx ' (64) 



o 



An for T = 



A 2 (T)=<( 2 (65) 
A N J . exp f 2 ^\ for T»T A 



12 T A "V. t aJ 

Here we defined A^ = 47re 2 n e /m and the Andreev temperature Ta = v-p/{2~Kd^). The extremely 
nonlocal form of the current-field relation (]64|) leads to an overscreening for A(T) < d^. This means, 
the magnetic induction B(x) changes its sign inside the normal metal, reaching for An *C c?n at the 
interface to the superconductor a value as large as —if/2, which points in opposite direction of the 
external field. The temperature is given by the Andreev temperature Ta, which plays the same role 
as the Thouless energy in the diffusive case. 
For the susceptibility we find 

X = ~474+12A 2 (T)/d N ' (66) 

In the limit A(T) <C (fo the susceptibility is 3/4 of — 1/47T, thus the screening is not perfectly 
diamagnetic. For A(T) ^> d^ screening is exponentially suppressed. Both lengths coincide, A(T) = 
cZn, at a crossover temperature T ~ 2Ta log(dj<s / X^s) , which can be considerably larger than Ta itself. 
Interestingly, it is exactly this temperature below which the first order transition in the nonlinear 
response occurs j74J]. The result for the linear susceptibility is displayed in Fig. || as the thin solid 
curve. 

3.2.3. Dirty limit 

Next we assume that the mean free path, I, is the smallest length scale (except for the Fermi wave 
length). In the absence of fields the Green's functions vary on a scale £n = (D/2ttT) 1 ^ 2 , and the 
dirty limit is realized if I <C £n- However in the situation with a space-dependent field, as seen from 
the kernel (^9|) , a stronger requirement is needed. The dirty-limit form of the current-field relation 
( pi]) is found only if vector potential changes slowly on the scale of the mean free path. Assuming 
this to be the case we obtain the susceptibility shown in Fig. ^| as the thin dotted curve. The typical 
energy scale, below which the screening saturates is given by the Thouless energy f/rh = D/d^ in 
analogy to the Andreev temperature Ta in the clean limit. 

In order to check whether the requirement of a slowly varying vector potential is satisfied, we 
define a local penetration depth X(x, T) as 

^y> 2 ( *). (67) 



In order to evaluate it we can proceed with the approximate form F w ~ exp[— x(2oj/D ) 1 ^ 2 ]. As a 
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Fig. 6. The susceptibility of the normal-metal layer for An = 0.003dN- The clean limit is indicated by a thin line 
reaching 0.75 for T — > 0. Shorter mean free paths , Z, (even as large as 10 4 cZn) lead to an enhanced screening at low 
temperatures. For shorter mean free paths, in the diffusive regime, a completely different temperature dependence is 
found. The dirty limit theory including a local current response deviates strongly from the correct result at higher 
temperatures. Inset: Comparison of experiment and theory. The thickness of the normal metal in the experiment was 
3.1/^m, corresponding to Ta = 540mK. The London length An = 22nm was taken from the literature. To fit these 
samples a mean free path I = 0.55cZn = 1.7/im had to be used. This is in rough agreement with the result I = 4/mi of 
a transport measurement. 



result we find for the local penetration depth 



Xnx/l if £ N (r)»d 



\K&(T)e°/t*F) /I if £ N (T)«d 



N 



For the local relation to be valid we need I < X(x,T) in the spatial region in which the screening 
occurs. For T <§C -&rh this means I -C X(d^), leading to the requirement I 2 <C An^n- The penetration 
depth A(cZn) is then given by X^d^/l. It is interesting to note, that this means that full screening 
A(cZn) <C g?n is only achieved for X^. For T ^> Ej-h screening takes place at x sa £n and we have 
I 2 <C An£n(T)- The low and high temperature conditions for the local approximation are different. 
Consequently it is possible that the dirty limit theory can be applied at low temperatures, whereas 
it fails to describe the correct behavior at higher temperatures. 

A comparison of the dirty-limit result and the more general result for a fairly short mean free 
path I = O.ld^, taking into account the nonlocal screening, is shown in Fig. ||. We clearly observe 
a difference between the two results for higher temperatures. From the above discussion it is clear 
that the temperature scale at which the deviation sets in is related to the Thouless energy i?xh, as 
demonstrated by the result in Fig. o. 
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3.2.4. Arbitrary impurity concentration 

Finally we allow for arbitrary values of the mean free path. A qualitative understanding may be 
gained from looking at the current vector potential relation in the limit I 3> c/n ■ hi the limit T <C Ta 
the zeroth-order Green's functions are given by the clean limit expressions. We approximate the 
kernel © by 

Since I 3> c?n, the exponentials may be expanded to first order. As a result, we obtain two contribu- 
tions to the current j(x) = j c iean + jimp(x). The first contribution jciean is given by the clean limit 
expression Eq. ( |64| ) . The second is an impurity induced contribution 

• r^ - 1 /"* N \x-x'\ + 2d n -x-x' ., ,., , 

= s^m^l — — A{x)dx - (70) 

This relations demonstrates when deviations from the clean limit become important. It is clear 
that the impurities cannot be neglected, if jimp (a;) is comparable to jciean- We estimate this by 
calculating the two contributions to the current using the clean-limit vector potential. Comparing 
the two contributions, we find that impurities can be neglected, if 

A c 3 ff (T)EA 2 (T);»4- (71) 

This equation defines a new length scale, the effective penetration depth A e ff. For the clean limit 
to be valid at T = the condition A c ff(0) > da has to be fulfilled, since in this case the screening 
takes place on the geometrical scale du- In the case A e ff(0) <C c?n the field is screened on a scale 
A e ff and the susceptibility is strongly enhanced in comparison to the clean limit. Nevertheless, the 
clean-limit behavior may reappear at higher temperatures, since A e g(T) grows with temperature. 

For T 3> Ta the deviations from the clean limit can be calculated as perturbation in the impurity 
scattering. The correction to g° leads to a finite superfluid density close to the superconductor via 
the factor 1 — g%{x') ss (£r/2i) exp(— 2x/£r) in the kernel. The range of the kernel is modified by 
the correction to /t°, which is non-negligible in this case. Since fl° ~ (£r/2Z) exp(— x/£t) we find 
for T>T A 

mUv x ,x,x>) = exp (l £ dx> ' « exp (2^) ■ (72) 

In deriving this equation the Dx-dependence was neglected by putting v x = vf- The range of the 
kernel is now given by £t, which is strongly temperature-dependent. For the current we find 

M ™ / dx ' e ~ 2(X +lX ~ X l)/iTA ( x ') » (73) 

again showing the importance of the length scale A e ff. In the limit A e ff(0) 3> £t the field cannot be 
screened on the scale ^, leading to a vanishing susceptibility. However, if A c ff (0) <C £t the field can 
be screened on a length scale smaller than £p and the susceptibility will be finite. 

For arbitrary mean free paths we have solved the screening problem numerically. The susceptibility 
as a function of temperature is shown in Fig. ^. The different curves refer to the clean limit and 
to mean free paths l/dys = 10 4 , 10, 1,0.1. In this plot we see the all the features discussed above. 
For large mean free paths the screening is enhanced in comparison with the clean limit at low and 
high temperatures, whereas the clean limit is restored in an intermediate temperature range. For 
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Fig. 7. A proximity wire as used in the experiment |l3| . The thickness is assumed to be <C £o and d is the 
distance between the strips 

small mean free paths (7 = O.IcZn m the graph) screening is suppressed also in the intermediate 
temperature range. 

A comparison to experimental data from Ref. |p"l|l is shown in the inset. Using the mean free 
path as the only fit parameter they can be reproduced well by the theory. Furthermore the fitted 
mean free path is in rough agreement with an independent experimental estimate obtained from a 
resistance measurement of the same sample above the critical temperature. 



3.3. Equilibrium supercurrent in a diffusive SNS junction 

The supercurrent through an SNS junction had been studied in the past mostly on the level of 
the Ginzburg-Landau theory [77| , which is valid near T c when the order parameter is small. Recent 
experiments entered the mesoscopic regime at low temperatures, where, e.g. the coherence length 
£n = y/D/2-KT is comparable to the sample size d. Here a more general approach is needed. 

Strong deviations from the Ginzburg-Landau behavior have shown in the current through a very 
thin (of thickness eN much smaller than £o)i diffusive normal wire on top of which an array of 
superconducting strips had been deposited E2[, see Fig. [?]. If there is a good metallic contact 
between the superconductors and the wire, we can assume the Green's functions under the strips to 
coincide f7§| with the values in the strips. These, in turn, can be assumed to be reservoirs, so the 
Green's functions are given by their bulk values The structure is thus equivalent to a chain of 
SNS junctions without barriers between the normal metal and the superconductor. 

We are interested in the equilibrium supercurrent. Hence, it is sufficient to use the Matsubara 
technique in the parameterization ([f9|). We thus have to solve the Usadel equation for the Green's 
function with the boundary conditions ( fi2| ) assuming <f>o = ±</>/2 at x — ztd/2. The general structure 
of the solution is such that 9 decays over a scale of £ w = y/V/2uj from the boundaries. If £ w <C d, 
the overlap between the induced 9 from the two superconductors is exponentially small. Then we 
can approximate the pair amplitude by a superposition of the solutions describing a semi-infinite 
normal part in the absence of phase gradients, i.e. 

F u (re) = e 1 ^ 2 sin 9+ (x) + e~ lv ' 2 sin 9~ (x) , 
tan(0±/4) - expl ±xd + A% I . (74) 




For T ^> E^h — D/d 2 this approximate solution is sufficient to calculate the current. From the sum 
( p^ ) we find I — I c sin 0, where the critical current is 

_ 64ttT x ^ d A 2 ^ f_d\ 
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Fig. 8. Critical current of the proximity wire compared 
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Fig. 9. The critical current at low temperature for 
different values of the gap, normalized to the Thouless 
Energy. A 2> &rh corresponds to long and A &rh 
to short junctions. 



At high temperatures the sum over Matsubara frequencies is dominated by the first term. If, fur- 
thermore, T C A, we arrive at EE, ' 



79 



64tt T d ( d 



3 + 2\/2ei? N £ N V 6 



The exponent reflects the fact, that the weakest part of the structure, i.e. the overlap region in the 
middle, defines the bottleneck for the supercurrent. 

As becomes apparent from the semi-logarithmic plot in Fig. || the critical current is approxi- 
mately described by the numerical fit, I c cx e~ T / T , with T* — YlE^ji:. This result arises because 
the convex exponential factor in (f76|), typical for diffusive metals, competes with the concave pre- 
exponent T 3 / 2 . This surprising result resolves a puzzle with the interpretation of the experiments of 
Ref. |D| . The observed linear T-dependence of the exponent is what one would expect in a ballistic 
system j77], |o), whereas the experiments were performed with diffusive metals. The present result 
agrees well with the experimental observation 

At lower temperatures, all Matsubara frequencies contribute to ( |T4| ) and we have to solve the 
Usadel equations numerically. The results are shown in Fig. |. As T < -&rh, the critical current 
saturates at a finite value I C (T — 0). The numerical results at very low T (Fig. |^ ) indicate that for 
a long junction, d 3> £o, the spatial decay of F limits the current to I c Rn oc £/rh. In contrast, for 



a short junction the gap is the relevant cutoff and I c Rn oc A |81|. This agrees qualitatively with 
estimates derived from the low-energy expansion of the Usadel equation J79[ The numerical results 
for the critical current under experimental conditions 113] explain the experiment quantitatively well 
(see Fig. |). 



4. Applications to nonequilibrium transport problems 
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4.1. Supercurrent under nonequilibrium conditions: The SNS Transistor 

In this section we will describe nonequilibrium transport problems. For this purpose the Keldysh 
formalism outlined before is particularly appropriate. As a first example we investigate how the 
supercurrent through an SNS structure is influenced by a nonequilibrium current in the normal 
part Jl5| . Specifically we study the set-up which recently has been investigated experimentally by 
the Groningen group |jl6^[l8"l| . The device is shown in Fig. [l^. Its center part is a normal conductor 
connecting massive reservoirs at voltages ±V/2. The distribution function in this part is obtained 
as a solution of the kinetic equation, which for a diffusive, mesoscopic normal metal reads 

Vdlf = . (77) 

Its solution 

f(E, V) = Q - |) r (E + eV/2) + Q + |) /° q ( E ~ el/ / 2 ) ( 78 ) 

has two temperature-rounded steps at the electrochemical potentials of both reservoirs. The step 
heights depend on the position along the wire; in this way the distribution function interpolates lin- 
early between the boundary conditions at y = ±L/2. This functional dependence had been verified 
in the experiments of Pothier et al. (l4|]. At the center of the normal conductor superconducting 
electrodes are attached, through which a supercurrent is flowing. By choosing a symmetric geometry 
we can assure that the effective chemical potential of the distribution function at the site of the 
superconducting contacts is zero. This guarantees that no net current is flowing between the nor- 
mal reservoirs and the superconducting leads. On the other hand, the nonequilibrium distribution 
function influences the properties of the SNS contact, and can be used to tune the supercurrent, 
which creates the possibility to use the device as a transistor |l6| . At low temperatures even a sign 
reversal of the critical current, i.e. a transition to a so-called 7r-junction has been predicted [fl5| and 
observed SQ. 




Fig. 10. Different realizations of the SNS transistors. The supercurrent is tuned by (a) a perpendicular or (b) a 
parallel normal current flow. 



To describe the system we have to evaluate the spectral functions as well as the distribution 
function. The former can be obtained from Eq. (pQj). The low-energy (E <§; i?Th) solution, calculated 
perturbatively, reads 

9 ~ -ztt/2 + (E/E Th )a(<P) + 0(E 3 ) 
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<j>x/d - (E/E Th ) 2 b{c/>) + 0{E 3 ) 



(79) 



with real- valued functions a and b. This indicates that the spectral supercurrent Imjjs}, vanishes 
in this energy range. At higher energies the problem can be analyzed numerically. However, more 
insight into the problem is obtained from an approximate analytic solution. For the calculation of the 
spectral functions between the superconducting contacts we neglect the widening of the geometry 
in the normal metal. In this case we can use the solution (|74|) obtained in the previous section. 
The approximation introduces errors of order one in numerical coefficients, equivalent to a change 
in an effective area. The result for the spectral supercurrent Im-jj^}, shown in Fig. O, displays a 



proximity-induced mini-gap 



of size E g ~ 3.2i?Th at = 0. This gap decreases with increasing 



i, when the induced pair-amplitudes from both sides start to interfere destructively, and vanishes 



at <f) = 7T | S3 1 _ At energies directly above the gap, Tmjjs} increases sharply, but rapidly decreases at 
higher E. At large energies, it changes sign and eventually oscillates around zero with exponentially 
decaying amplitude. 
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Fig. 11. The spectral current Im(jE) as a function of 
energy for different values of the phase difference <j>. At 
higher energies we find oscillations as expected from eq. 



Fig. 12. The supercurrent as function of control volt- 
age (at T = 0) and temperature (at V = 0) for various 
values of rf>. 



To evaluate the physical current (42) we combine the spectral supercurrent and the quasiparticle 
distribution function / between the reservoirs at electrochemical potentials ±eV/2 derived above. 
Although no net current is flowing between the normal metal and the superconducting electrodes, 
the coupling might influence the energy dependence of the distribution. We will now show that this 
is not the case: Compared to the distance between the reservoirs L the proximity effect is confined 
to a relatively narrow region ds- For a quantitative analysis we have to distinguish between the two 
components /l/t> governed by the kinetic equations and (|35|). For the supercurrent we 
need to determine only /l. Because of the symmetry of the setup, the component /t vanishes at 
the site of the superconducting leads. 

As the voltage is applied between the reservoirs, the gradients of f L / T point in y-direction, i.e. 
perpendicular to je which points in x-direction. As a consequence, the scalar products V/l/t • ]e 
vanish and (|34j) and (|35|) decouple. Moreover, /l satisfies the same boundary conditions (|3^) at 
the two reservoirs at ±V/2. Altogether, (p3) implies that /l has a constant value given by (|3S 
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throughout the wire. It is easy to check, using 2/ = 1 — /l — /t with fx = and (p3|), that the 
total distribution combines to the double-step form described before. 

At first glance, it may look surprising that /l has such a strong nonequilibrium form in the normal 
metal next to the two superconductors, which are assumed to be in equilibrium |56fl . It turns out to 
be a consequence of the symmetry properties of Andreev reflection. As the SN-interfaces are good 
metallic contacts, the boundary conditions eq. p5| ) guarantee, that at the interface 2?l coincides 
with the value in the superconductor. Thus at \E\ < A we find from eq. ( |39|) that 2?l = 0. This 
implies, that the SN-interfaces decouple /l in the normal metal from its value in the superconductor. 
A well known physical consequence of this property is, that below the gap, when quasiparticles are 
transmitted by the Andreev reflection, no heat is transferred || . The operation of the SNS-transistor 
without tunneling barriers is based on this very idea [fl6| . 

Our treatment of the 2D system is approximate. A more complete solution would also take the 
bending of supercurrcnt lines into account, which leads to an effective renormalization of the cross 
section area. This problem has been solved by Volkov [^6| for a related system with small ds and 
tunneling barriers at the SN-interface. An analysis of his solutions confirms the aforementioned 
argument, that the effective area of the normal region between the superconductors is larger than 
djs[d. This renormalization is reduces at higher temperature and vanishes at T> £?Th- 

The form of the distribution function implies that in the energy window — V/2 < E < V/2 current- 
carrying-states are depopulated and the supercurrent ( f42| ) is strongly suppressed. For V < E g this 
has no effect because of the gap. At higher voltage the supercurrent decays rapidly with increasing 
voltage (cf. Fig. |l2| ). This is exactly what has been observed in the experiments |l6|]. Furthermore, 
at still larger voltage eV > 10i?Th the dominant contribution to the integral in (|42|) outside this 
energy window is negative (see Fig. [ll 12), so the total supercurrent changes sign. We thus find a 



transition to a 7r-junction |p4| , controlled by nonequilibrium effects. This prediction has been verified 
in recent experiments |l7],|l8|]. The effect is rather pronounced, the critical current of the 7r-junction 
is approximately 30% of I c at T = eV = 0. 

For high voltages or temperatures eV,T ^ E s we find a closed expression for the current (p2[). 
The integral is dominated by the poles of /l, given by (|33|), at imaginary frequencies E = iuj n ±V/2. 
Also je can be found by analytic continuation of the approximation (j7J) . We thus find Is — I c sin <f>, 
with critical current 



IcRd = _^ e -V^: C0S ( ; v ( , 

3 + 2v^ \2y/n v ETh J 1 T ( } J e ) forT»£ g 




for eV> E K , T < E g 

(80) 



Here n v = ttT + ^{irT) 2 + (eV/2) 2 and 

f vr/2 for eV > E z , T < E g 

n for T»E g ■ ( 81 ) 

This result indicates that with increasing temperature the onset of the 7r-junction behavior is weak- 
ened and shifted to higher voltages. Our numerical results, shown in Fig. [H^, demonstrate that it is 
still pronounced as long as T < £rh- 

Above we have discussed the limit of A > Exh ■ For smaller values of A the minigap E g is limited 
by A rather than E^h, see Fig. [l4| However, the qualitative behavior, including the transition to 
the 7r-junction regime persists. 




Fig. 13. The supercurrent as a function of the control 
voltage for various values of T at <p = tt/2. 



Fig. 14. The supercurrent as a function of the control 
voltage (at T = 0) for various values of A/&ph at </> = 
tt/2. 



4.2. Dissipative current in a mesoscopic wire 

Perhaps the simplest system displaying the influence of the proximity effect on transport properties 
is shown in Fig. IB]: A mesoscopic wire connects a superconducting and a normal reservoir. At the 
interfaces additional barriers (indicated by r in Fig. |l5|) may exist. In this section we assume no such 
barriers to be present. We assume the superconductor is assumed to be at voltage V = 0. As the 
dissipative current is a genuine nonequilibrium quantity we again will use the Keldysh technique. In 
addition to the previously used boundary condition at the superconductor, we impose = at the 
normal reservoir. Hence, the spectral supercurrent je vanishes at the normal lead. As the retarded 
Usadel equation (|5(]) conserves this spectral current, this implies that je = everywhere and hence 
d x \ = 0. This means that no 'supercurrent' flows in this system. On the other hand, the 'normal' 
current j n is influenced by superconducting correlations, and we have to study their behavior in an 
electric field. 

Approximate solutions of (|5^) can be obtained by the same methods as described above. We will 
now concentrate on the case £>rh <C A. For E <C i?Th, iterating ( |50| ) from the E — 0-solution yields 

8 E r fX"K\ , 7T 



7T 2 E Th ' 



in( T )]-^ (82) 




Fig. 15. A thin normal wire between normal and superconducting reservoirs. 
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where x = x/d. For E 3> _Exh we have 



tanh I ~ I = - tanh ( ^ ) exp[fc(a; - d)] (83) 



4 I V 8 J 



with k = yf— 2iE/T> . For the calculation of the dissipative current we only need the distribution 
function /t- As je = 0, the kinetic equations decouple and (Q) reduces to a simple spectral diffusion 
equation. Relating /t to the potential and T>i to the conductivity, we note that our expressions ( |43| ) 
and (|34|) are analogous to Ohm's law and the continuity equation on a spectral level, respectively. 
It is now straightforward to integrate ([mJ) and to calculate the current (fl3|), from which we obtain 
the differential conductance 

/°° dE 
5 9(E) ■ (84) 
-co 2Tcosh 2 [(£;-i/)/2r] yv ; 1 ' 

Here the spectral conductance g(E) is defined by 

5(£;)=GNd (r^)) • (85) 

Again this expression is the spectral form of Ohm's law. Using the approximate solutions (|82|) and 
(S3), we calculate 2?t and g(E), and arrive, e.g. at V = 0, at 

g _ 1 | h T 2 / E d for T«E Th 

G N UVVT forT»^ Th ^ j 

with a w 0.049 and 6 w 0.42. Results are displayed in Fig. |l^. The conductance approaches the 
normal-state value for T 3> E^h, because in this limit the proximity effect is confined to a short 
portion of the wire of length £n- However, also at T — > the correction to the conductance vanishes. 
This has been first noted by Artemenko et al. |85) for short contacts, and later in longer wires by 



Nazarov and Stoof |19j,|2C|, 86 1. The theoretical prediction was confirmed in experiments |2l|,|2^]. To 
the best of our knowledge, no simple physical explanation of this surprising coincidence has been 
provided yet. The fact that the conductance is unchanged does not imply that the proximity effect 
is absent. Even at E = the solution yields F R ^ 0, and the density of states Re(G R ) displays 



a proximity-induced gap-structure, as shown in Fig. 17. 

At E = 0, the local spectral diffusion coefficient T>t coincides with the normal-state value every- 
where, see Eq. (|S2|). Returning to the definition ( |37j ) and using the parameterization (^9| ) we find 
here @ 

X> T = (ReG R ) 2 + (ImF R ) 2 = cosh 2 (Re(6>)) > 1 . (87) 

This means, that it is a combination of the density of states and the Maki-Thompson paraconduc- 
tivity term [[f7|, |S7f . Both terms compete. At finite E, the Maki-Thompson contribution is stronger 
than the suppression of the DOS JT^ . At E — both contributions, although each of them varies 
in space, add up to 1. 

In order to study the space-dependence more thoroughly, we look at the local DOS in our model. 
Our results (|l7]) indicate, that due to the induced superconductivity, a pseudo-gap in the DOS 
opens, i.e. the DOS is in general suppressed but does not vanish at low energies. If we move away 
from the superconductor, the pseudo-gap is weakened and the DOS approaches its normal-state 
form when we are close to the normal reservoir. 
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Fig. 16. Linear conductance of a proximity wire of 
length d with transparent interfaces. The energy is mea- 
sured in units of the Thouless energy &rh = D/d 2 as- 
sociated with the length d of the wire. 



Fig. 17. Local density of states in a proximity wire 
with transparent interfaces for A ^> &rh a t distances 
x/d = 0, 0.1, . . . , 1.0 (top to bottom at low E) from the 
normal reservoir. 



4.3. Supercurrent and interference effects on the dissipative current 



We will extend our study of the dissipative current to a situation where also a supercurrent can 
flow and therefore interference of the induced correlations starts to play a role. This type of effects 
was first studied experimentally p5[ | in cross-shaped diffusive metal structures with different super- 
conducting leads closed by a SQUID loop, later also in analogous semiconducting structures p8[ , 
as well as in metallic structures with an Aharonov-Bohm loop [^,^9|. Usually, these structures are 
referred to as "Andreev interferometers" . The common key observation is, that the conductance 
depends periodically on the magnetic flux penetrating the loop with a period of a superconducting 
flux quantum O = h/2e. Moreover, it turned out that these oscillations could be detected even if 
£n <C d where d is a typical geometrical sample dimension. I.e., the the proximity effect had a longer 
range than expected from simpler arguments. 

These effects have been extensively studied theoretically ]l9|,|^,^l],|6,90,9^]. We want to outline 
the main concepts in these systems, considering as an example the experiments of Ref. [p3[ . The 
layout of the system is shown in Fig. fL8l Applying the arguments of the previous chapter, we 
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Fig. 18. A proximity wire containing a loop. The cross-section areas of the respective components are Ai 
the field-induced supercurrent in the ring part. 
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recognize that the supercurrent is confined to the loop. Instead of dealing with the vector potential 
in the derivatives in (pfj), we remove it by a gauge transformation \ = xo — 2e J Q R dr' A(r') where 
the integral runs around the loop. This implies, that the vector potential produces an artificial phase 
difference of (j> = $/§o (modulo 2ir) at one point in the loop, which we choose to be point B. This 
procedure already guarantees the /i/2e-periodicity of the results. 

We know from the previous sections the form of the Usadel equations in the single branches. In 
the present problem we also have to account for the nodes. To do so we apply the idea of Zaitsev |)2| : 
By regarding the Green's functions in the connecting wires as cross-sectional averages of the full 3D 
Green's function, which have well-defined derivatives in space, we will conserve the "spectral flow" . 
This means 

Gn Gin 

for all branches n and m and at nodal points we have 

d 



G- — G = 

OX n 



(89) 



where n labels the branches, and the derivatives point in the direction of branch n away from 
the node. This approach has been used by Nazarov |)3,94| to derive a simple circuit theory for 
T = V = 0. Here, we will assume that the branches of the loop have half of the cross-section area 
of the other branches. 

Only very few analytical results have been derived for this system |l^,|9^]. Here we will focus on 
numerical results. As shown in Fig. [l9[ the spectral conductance shows a non-monotonic shape as 
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Fig. 19. Flux-dependent spectral conductance of the 
loop structure, assuming A 1/4 = 2A 2 /3 all di = d and 
A> E d = D/d 2 . 



Fig. 20. Linear-response conductance of the loop as 
described in the caption of Fig. Il9l 



in the previous section. Provided that the induced correlations persist from the superconducting 
electrode to the loop, the low energy results, for E <C i?Th, depend on the flux, whereas the 
conductance at E E^h is flux-independent, see Fig. [l9| This agrees with the intuition on the 
range of the proximity effect developed in the previous sections. In the loop, the vector potential 
causes the correlations to interfere destructively, so the induced conductance correction is suppressed 
by the presence of the loop. At E 3> E^h, the correlations have already decayed before they reach 
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Fig. 21. Full period of normalized conductance oscillations in the loop structure at V = 0. 




the loop and are therefore not influenced by the flux. On the other hand the temperature dependent 
conductance also depends notably on the flux at T > £rh, see Fig. so we indeed recover 
the observed long range proximity effect there [^3-25|. These seemingly paradoxal observation can 
be resolved by examining (|84|). In order to calculate the conductance Fig. ^o] from the spectral 
conductance Fig. [l9| we have to convolute the latter with the derivative of the Fermi function. This 
distribution function is peaked around E = and has the width T, i.e. the window of size -E/rh 
around the Fermi surface where the spectral conductance is actually flux-dependent with a relative 
weight of -Bxh/7 1 . Thus, the conductance depends sensitively on the contribution of this low-energy 
range, where the proximity effect has a long range p3fl , which is only limited by the phase-breaking 
length 1$. 

We note that the intuition on £n developed by regarding thermodynamic quantities as the su- 
percurrent in an SNS junction (section |3.3| ) fails when describing nonequilibrium quantities. On the 
other hand, we have shown that the same microscopic mechanism (characterized by the same Green's 
functions) may lead to the exponential temperature dependence of equilibrium quantities as well as 
to the "long range" tail with power-law temperature dependence of nonequilibrium phenomena. 



4.4. Transport through tunneling barriers 

So far, we only considered situations, where the current flow does not have to overcome tunneling 
barriers. We found, that the local spectral diffusion coefficient T>t = cosh 2 (Re6>) > 1 and thereby the 
global conductance is always enhanced in comparison to the normal state. However, as we explained 
before (section |3.l| ), the tunneling conductance is proportional to the density of states, i.e. at low 
energies, E < E^h, it is suppressed in comparison to the normal state as a manifestation of the 
induced superconducting gap In general, these two mechanisms compete. Moreover, due 

to the direct coupling to the normal reservoir, correlations immediately feel the presence of normal 
electrons and only a pseudo-gap develops, which turns into a real minigap when the reservoir is 
decoupled from the system via a tunneling barrier [|l9[ ||| . 

Following the lines of Ref. ]96| with boundary condition (|54|), we can calculate the distribution 
functions and from there the current and conductance. The result has the form (|84|), however, the 




Fig. 22. Linear conductance as a function of temper- Fig. J23. Linear conductance of the interferometer in 

ature at different tunneling resistances within the diffu- Fig. E3 as a function of the phase in the tunneling- 

sion/tunneling crossover at <f> = 0. At (j> = ir we have dominated regime; inset: diffusion-dominated regime. 

G = Gn. The inset shows the structure used for cal- 
culations. All arms 1-4 are assumed to have the same 
lengths d and cross sections A. 



spectral conductance now accounts for the series of the tunneling resistance and the diffusive metal 

Here rx = Rt/Rn, where i?N is the normal state resistance of the wire, and Vt — Re{G R } is the 
reduced DOS at the interface. This result implies that already a rather low barrier of size ~ i?N 
gives rise to a tunneling-like structure. In typical experiments, this is of order 1 . . . 10O. In Fig. [2^ 
we show the conductance of a structure relevant for experiments p5[ demonstrating that at low T 
for tt > 1 the conductance is reduced compared to the normal state value. 



As we showed in section 4.3, a magnetic flux through an Andreev interferometer causes interfer- 



ence effects which influence the conductance. This interference turns out to be destructive. For the 



system without barriers studied in section 4.3, this manifests itself in the fact that the conductance 
oscillations start from a zero-flux maximum. If the interferometer contains tunneling barriers, we 
expect that the interference effects inside the interferometer are qualitatively not affected, however 
the conductance eq. ( |90| ) is dominated by tunneling into the induced gap. The destructive inter- 
ference will weaken induced superconductivity, therefore lift the gap and consequently increase the 
conductan ce. Thus in this system we also expect magneto-conductance oscillations as described 



in section LS , however, they will depa rt f rom a zero-flux minimum. By combining the theory of 
the interferometer developed in section [4.3| with the conductance formula (|9^) , we can study these 
structures and confirm this picture |95|, |97|. In Fig. ^3] we show a half-period of the conductance 
oscillations at different temperatures for a tunneling-dominated system in comparison to the same 
system probed through metallic contacts which have a different sign of the flux-dependence. This 
effect has recently also been observed in experiments |9q|. 
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Fig. 24. A quasi-2D proximity film. The contacts 
A and B are used as voltage and C and D as cur- 
rent probes. An alternative setup: A and C are voltage 
probes, while the current flows through B and D. 
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Fig. 25. A qualitatively equivalent circuit with the 
probe configuration as in Fig. E4. 



4.5. Four-Point conductance 

In the previous sections we argued, that the proximity effect increases the conductivity of a normal 
metal wire compared to the normal state. If the current is probed through tunneling barriers, this 
increase competes with the opening of a superconducting gap which can cause the total conductance 
to decrease. However, a decrease in the conductance was also observed in experiments (24],|2[| where 
the current does not have to pass through tunneling barriers. This effect cannot be explained by 



the non-monotonic conductance described in section 4.2, because despite its T-dependence, the 
correction to the normal state conductance is still non-negative. 

In the experiments |pi],p5||, the current was probed by a particular four-probe configuration on a 
sample of appreciable width. Hence the current flow cannot be described by simple quasi-lD models 
as above. We will now develop a quasiclassical description of such measurements and show how the 
inhomogenous conductivity of the sample leads to a decrease of the four-point conductance. 

Let us consider a normal film connected to a superconductor with four point- like contacts, which 



lead into measuring reservoirs after a distance much larger than l^, see Fig. 24. In this configuration, 
the probes do not influence the spectral properties, so the solution of the retarded Usadel equation is 
analogous to the quasi-one dimensional case with an infinite tunneling barrier separating the normal 
reservoir from the sample 

d x 9(x = d) = 0. (91) 
At large energies, 9 is small and we reproduce (|83|). At low energies, a perturbative expansion gives 

e E (x) = ^j(2-j)-tn/2. (92) 

The kinetic equation is now two-dimensional where the current contacts are represented by source 
terms 

V(2>r(r)V/r) = H^r - r c ) - S(r - r D )) , (93) 
A "no current flow" condition at the normal metal edges yields 

d n h = ■ (94) 

A formal solution of Eq. (|34|) reads 

f T (E,r) = I E (G E (r,r c ) - G E (r,r D )) , (95) 
where Q E = (X> T V + X> T V 2 ) -1 is the Green's function of the operator (||). Making use of |33| ||) 
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Fig. 26. The linear conductance G(T) for films of different widths w/d = 0.05,0.5, 1.0 calculated for d = L and 
T c = 5.7 10 5 &rh- Left inset: The same curve at w = 0.5d. The T-axis is zoomed to demonstrate the presence of a usual 
ID-type non-monotonic behavior at T ~ i?Th- Right inset: G(T) for a wide film and T c = 50&rh- The amplitude of 
the negative conductance peak is increased due to the effect of a superconducting gap A(T). 



and (|95|), and integrating Ie over energy we obtain the total current / and arrive at the expression 
for the differential four-point-conductance G = dl/dV: 

where 

r.BC r, BD rAc \ r. AD 

g{E) - G ^gBC^gB D _gAC + gA D (97) 

is the spectral conductance. We introduced the notation = G(ri,rj) and Qo is the Green's function 
of ( [34]) in the normal state {Me{t) = <th). The spectral conductance (|97]) calculated numerically 
from Eqs. @, (|9|) and (||) is presented in Fig. 

For narrow films the well-known results of quasi-lD calculations Jl9|,|20| are qualitatively repro- 
duced: the linear conductance G(T) exceeds Gn at all temperatures, showing a non-monotonic 
feature at T < (for simplicity we put L = d here and below). sFor broader films G{T) decreases 
below the normal-state value at high temperatures and reaches the minimum at T ~ 10i?Th- At 
lower T the conductance grows with decreasing T, becomes larger than Gn and then decreases again 



down to G(T = 0) = Gn similarly to the ID case (see the left inset in Fig. 26). The behavior of 
g{E) = G(E, T = 0) as a function of energy (voltage) is qualitatively identical to that of G(T), the 
negative peak at E ~ lOi/rh turns out to be even somewhat deeper. 

In order to develop a semi-quantitative understanding of these results, let us consider a wire 
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network representing the interconnections between the probes as in fig. ^5|. A similar equivalent- 
circuit model was previously used for a qualitative description of inhomogeneous superconducting 
films |Q. Exploiting the analogy between /t and the electrical potential in a conventional circuit 
(see above and p6[]), Kirchhoff's laws for the spectral conductances can be derived Jl9| , |92| ]. For the 
present circuit Sig. ^5|, we find (c.f. [p8[) 

4 

SNct = 9394^9^ , (98) 

i=l 

where the gi are the spectral conductances [^9[Q of the wires 1-4 calculated analogous to previous 
sections. At T > i/ph only the wire 1 directly attached to a superconductor acquires supercon- 
ducting properties, whereas the proximity effect in the wires 2, 3 and 4 is suppressed. Thus only g\ 
increases, and (72,3,4 remain unaffected. According to Eq. (|9|) g^ ct decreases below Gn- At T w i?Th 
the proximity induced superconducting correlation penetrates into all four wires, i.e. 52,3,4 increase, 
leading to the increase of (/Net above Gn- Close to A, the conductance correction has a peak cor- 
responding to the peak in the DOS. If close to T c , A becomes close to E^h, the conductance of 
interconnection 1 is peaked and the whole effect becomes even stronger, in agreement to our numer- 
ical result. 

A semi-quantitative analysis |2(| shows, that the crossover energy between excess and deficit con- 
ductance is of the order of max(D/d 2 , D/W 2 ). A careful analysis of the experiment s p4] , |25| shows 
that this explains the different signs of the conductance for different measurements f26| . 



5. Conclusion 

We outlined the concepts of quasiclassical Green's functions. They provide a unified view of 
ballistic and diffusive systems in equilibrium and nonequilibrium situations. We applied the theory 
to describe mesoscopic proximity structures. 
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B. J. van Wees. This work was supported by the DFG through Sonderforschungsbereich 195 and a 
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